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The current - voltage dependence of the transverse tunnehng current through the electron or hole 
traps in a DNA is investigated. The hopping of the charge between the sites of the trap and the 
charge-phonon coupling results in a staircase structure of the I-V curve. For typical parameters 
of the DNA molecule the energy characteristics of a DNA trap can be extracted from the I-V 
dependence, viz., for a small gate voltage the phonon frequency and for a large gate voltage the 
hopping integral can be found from the positions of the steps in the I-V curve. Formation of the 
polaronic state also results in the redistribution of the tunneling current between the different sites 
of the traps. 
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For the past several years intense research on electron 
(or hole) migration through DNA has established that 
these molecules are capable of transportingcharges over 
a distance of at least a few nanometers 

Mi a This 

is primarily achieved by oxidation of guanines (G) which 
generates a guanine radical cation. Guanine has the low- 
est oxidation potential of the common DNA bases 
and therefore a guanine radical cation can only oxidize 
another G. Stacked Gs such as GG, or a GGG, having 
much lower ionization potential than that of an isolated 
G, known to act as hole traps ^J. Since accumulation of 
holes at G-sequences usually leads to deleterious effects, 
including mutations charge migration through DNA 
may have important biological consequences. Electronic 
propreties of both electron and hole traps depend cru- 
cially on the charge dynamics inside the trapping spots. 
Therefore, it is very important to extract the parame- 
ters which characterize such dynamics. Here we show 
how the energy characteristics of a DNA trap can be ex- 
tracted from the current (I) - voltage (V) dependence of 
the transverse tunneling current through the electron or 
hole traps. Hopping of the charge between the sites of the 
trap and the charge-phonon coupling results in a stair- 
case structure of the I-V curve. For a small gate voltage 
the phonon frequency and for a large gate voltage the 
hopping integral can be found from the positions of the 
steps in the I-V curve. Formation of the polaronic state 
also results in the redistribution of the tunneling current 
between the different sites of the traps. 

Since our interest is in the local properties of DNA 
traps, transport in the direction perpendicular to the 
backbone axis (transverse transport jal) is important. In 
this case, if the electrodes have a relatively small width 
tunneling occurs through a single DNA base pair. The 
saturated (unstructured) tunneling current then depends 
on the particular type of base pair 's^ . This fact can be 
used to discover the sequence of DNA by scanning it with 
conducting probes. We demonstrate here that not only 
the saturated value of the tunneling current but also the 
structure of the I-V curves can provide important infor- 
mation on the properties of the DNA, in particular about 



the trapping spots. This is because the tunneling cur- 
rent through the system is determined by the density of 
states (DOS) of the system. For a finite system the DOS 
has peaks corresponding to discrete energy levels. These 
peaks will result in a staircase structure of the tunnel- 
ing current as a function of the applied voltage whenever 
the Fermi levels align with a new state of the system 
and thereby open an additional channel for tunneling. 
Therefore, from the staircase structure of TV curve one 
can learn about the energy spectra of the system. For 
DNA the trapping spots consist of a finite number of 
base pairs. Hopping between the base pairs within the 
traps determine the energy spectra of the spots. In ad- 
dition to the energy scale due to hopping there is also an 
energy scale due to the electron-phonon (or holc-phonon) 
ineraction. Finally, for DNA trapping spots the I-V de- 
pendence has two types of staircase structure, one due 
to hopping and the other due to the phonons. We have 
explored the interplay between these two effects. 

The tunneling transport through a single molecule or 
a quantum dot with electron-phonon coup ling has been 
extensively studied in the literature j^, llfl HJI ■ The 
main outcome of these works is the staircase structure 
of I-V curves due to phonon sidebands. The heights of 
the steps in this structure depend on the strength of 
the electron-phonon interactions, temperature, and on 
the equilibrium condition of the electron-phonon system. 
These studies were mainly restricted to a molecule with a 
single electron energy level, although a general approach 
to a many level system is also formulated . 

The DNA trap can be considered as a system of a few 
molecules (base pairs) with hopping between them and 
the electron-phonon coupling. Then in the I-V curves we 
should observe the interplay between the staircase struc- 
ture due to hopping between the molecules and due to 
the phonon sidebands. Since the tunneling occurs only 
through a single base pair the I-V structure should also 
depend on the position of the base pair through which 
the tunneling current is measured. Here we consider only 
the hole traps and the tunneling current of holes, but the 
analysis is also valid for electron traps and electron trans- 
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port. Whether it is a hole transport or electron transport 
depends on the gate potential, i.e. on the position of the 
chemical potential at zero source-drain voltage Vsd- 

We disregard below the effects related to a Coulomb 
blockade or to a double occupancy of the DNA traps, 
assuming that the repulsion between the holes is strong 
enough, although for some range of parameter of the 
traps the bipolaronic effect may become strong and a 
two-hole system can have lower energy than a single-hole 
system For a single hole in the trap the Hamil- 

tonian of the trap and the electrodes consists of three 
parts: (i) the trap Hamiltonian which include tight- 
binding hole part with hopping between the nearest base 
pairs, the Holstein's phonon Hamiltonian with diagonal 
hole-phonon interaction , (ii) the Hamiltonian of two 
leads, left (L) and right (R), and (iii) the Hamiltonian 
corresponding to tunneling between the leads and DNA 
traps 



trap 



'Hleads + 'Ht, 



(1) 



where 



-t- 



'Htrap = £ a\ai - ^ [al"i+i + h.c. 

i—l i 

(3) 
(4) 



'Hleads — ^ ^fe ^i,fc^a,fc, 
k,a=L,R 

a=L,R,k 



where is the annihilation operator of the hole on site 
(base pair) i, e is the on-site energy of the hole in the 
trap (same for all base pairs within the trap and is de- 
termined by the gate voltage or doping of DNA), bi is 
the annihilation operator of a phonon on site i, t is the 
hopping integral between the nearest base pairs, u is the 
phonon frequency, x is hole-phonon coupling constant, 
and da,k is the annihilation operator of a hole in the lead 
a — L, R with momentum k. The index z = 1, . . . , iVj in 
Eq. (0) labels the sites (base pairs) in the trap and Nt is 
their total number. Tunneling from the leads to the trap 
occurs only to the site iq with the tunneling amplitude Iq. 
In the hole-phonon part of the DNA Hamiltonian Titrap, 
we include only the optical phonons |l6j with diagonal 
hole-phonon interaction. 

We describe the process of tunneling through the trap 
as a sequential tunneling '17] . In the weak lead-trap cou- 
pling regime the tunneling Hamiltonian Tit can be con- 
sidered as a perturbation which introduce the transitions 
between the states of the trap Hamiltonian, Titrap- We 
denote the eigenstates of the trap Hamiltonian without 
couphng to the leads as |0,to) with energy E^o.m for the 
trap without any holes, and |l,n) with the energy 
for trap with a single hole. In the weak lead-trap cou- 
pling limit the master equation for the density matrix of 




(b) 










f 











FIG. 1: Current vs the sourse-drain voltage shown for two 
base pairs in the trap (Nt = 2) and different values of phonon 
frequency and hole-phonon interaction strength: (a) 7 = 0.1, 
A = 0.5; (b) 7 = 0.1, A = 1.0; (c) 7 = 0.2, A = 0.5; (d) 
7 = 0.2, A = 1.0. Black line corresponds to e = 1.3t while the 
red line is for e = 2.7. 



the trap reduces to the rate equation [T^ for probability 
Po,m to occupy the state |0, m) and probability Pi_„ to 
occupy the state |l,n), 



dt 



L,R r 



1 



(5) 



dP, 



0,m 



dt 
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7i,a—L,R 



Q.nm li'i 



E 

n,a—L,R 



0,m 



-^0,m E -^0'' 



(6) 



The distributions P^l^ and PqI^ are the corresponding 
equilibrium distributions with temperature T, P^l^ — 
exp(-Si,„/fcT)/^„,exp(-Pi,„VA:T) and P^^^ = 
exp{-Eo^m/kT) /J2m'^^Pi~^o,m'/kT). Here r is the 
relaxation time which is assumed to be same with and 
without a hole in the trap. The transition rate W^^^ is 
the rate of hole tunneling from the state |1, n) of the trap 
to the a — L,R lead leaving the trap in the state |0, m). 
Similarly, the rate W^~^„ is the rate of hole tunneling 
from the lead a to the state |l,n) of the trap, while orig- 
inally the trap was in the state |0,m). These rates can 
be found from Fermi's golden rule 

W^i:;™ -ro/„(i?i,„-i?o,,„)|(0,m|a,„ \l,n)\\ (7) 



w^Zi = To [1 - U {El,. 



Eo.m)] |(0,m|ai(, |l,n)| 



(8) 

where Fq = 2TTtop/h and p is the density of states in the 
leads, which is assumed to be the same in "L" and "R" 
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FIG. 2: Current vs the sourse-drain voltage shown for three 
base pairs in the trap {Nt — 3) and 7 = 0.2, but for different 
tunneling points io and different values of hole-phonon inter- 
action strength: (a) io = 1; A = 0.5 and e — 1.7t (black), 
e = 3.0t (red), A = and e = 1.7t (green); (b) io = 2: A = 0.5 
and e = 1.7t (black), e = 3.0t (red), and A = 0, e = 1.7t 
(green); (c) io = 1: A = 1 and e = 1.7t (black), e = 3.0t (red), 
(d) io = 2; A = 1 and e = 1.7t (black), e = 3.0t (red). 



leads, and fai^) is the Fermi function of the lead a with 
a chemical potential ^a- 

For the stationary case the time derivatives of Pi,„ 
and Po,m are zero and Eqs. (|SJ)-||HJ) becomes the sys- 
tem of linear equations with condition of normalization 
X^n-^i," + J2m^o,m = 1- Then the stationary current 
can be calculated as 



rn ^' L,mn 



Pi, 



I** L,nmj 



(9) 



We have calculated the current numerically as a func- 
tion of Vsd for different values of on-site energy, e, which 
can be changed by the gate voltage or by doping. By 
varying Vsd we are keeping the on-site energy e the same 
and vary the chemical potentials of leads as /^l = Vsd /2 
and hb. = -Vsd/2. 

To determine the tunneling current we first calculate 
the energy spectra of the DNA Hamiltonian [Eq. ^] 
without holes and with a single hole in the trap. The 
Htrap then is just the Hamiltonian of free phonons at 
each site of the trap. With a single hole we make the sys- 
tem finite by introducing limitations on the total number 
of phonons 19] Uph^i < 10, where Uph.i is the number 
of phonons on site i. The eigenf unctions and eigenvalues 
of Ti-trap can then be found numerically and the corre- 
sponding transition rates Eqs. Q-ljHJl can be calculated. 
Finally, we solve the system of linear eqiations ®-® 
for a given bias voltage, and substitute this solution into 
Eq. to find the stationary tunneling current. 

There are five dimensionless parameters which charac- 
terize the I-V structure: the nonadiabaticity parameter 
[l8| 7 = huj/t, with a typical value of ~ 0.01 — 0.5 for 
DNA, the canonical hole-phonon coupling constant ^3 
A = x^/ (2/ia;t) which is ^ 0.2 — 1 for DNA, dimensionless 
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base pair index, / 

FIG. 3: The average number of holes (a) and the average 
number of phonons (b) for a single hole system in a GGG trap 
are shown as a function of the base index. Dots and triangles 
corresponds to hole-phonon interaction strength A = 1.0 and 
0.5 eV respectively. 



bias voltage Vsd/t, on-site energy e/t, and the ratio of the 
relaxation time and the tunneling time tTq. 

The calculations have been performed for Nt — 2 and 
Nt — 3, i.e., for 2 and 3 base pairs in the trap. The 
example of such a system could be the guanine hole traps: 
GG and GGG spots surrounded by adenines. In all the 
calculations we kept the ratio of relaxation and tunneling 
time equal to 1 {tTq = 1), i.e., the hole-phonon system 
in the trap is not in the equilibrium. Different values of 
tTq, ranging from tTq <^ 1 (equlibrium case) to tTq ^ 
1 (unequlibrium case) do not modify qualitatively the 
behavior of the I-V curve. To detect the phonon steps 
in the I-V curve the temperature should be less then the 
phonon frequency and so we keep the temperature equal 
to O.Olt. 

In Fig. 1 our results are shown for two base pairs (sites) 
in the trap. The tunneling occurs through one of the 
site, io = 1. For uncoupled hole-phonon system the I- 
V dependence has two steps corresponding to two single 
hole energy levels. The distance between the steps is 
SVsd = 2t. For a small hole-phonon coupling constant 
A = 0.5 [Figs. 1 (a,c)] the additional structures of width 
SVsd — fi^ due to the phonon sidebands appear only at 
the first step and the second step can still be clearly dis- 
tinguished. At the same time for a large gate voltage 
(large on-site energy e), the phonon steps are suppressed 
and the I-V structure becomes similar to that of a zero- 
coupling strength, which is shown in Figs. 1 (a,c) by red 
lines. For a strong hole-phonon interaction (A = 1) the 
phonon steps suppress the steps due to inter-site hop- 
ping within the trap [Figs. 1 (b,d)]. This suppression 
becomes stonger for a larger non-parabolicity 7, which is 
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illustrated in Figs. 1 (b,d) by a black line for 7 = 0.1 and 
7 = 0.2. With increasing gate voltage the phonon steps 
disappear and the I-V curve shows a clear two-step struc- 
ture. The origin of such a suppression of phonon steps is 
the following. If the relaxation time is much smaller than 
the tunneling time then before tunneling the hole-phonon 
system will be in equilibrium. At low temperatures this 
means that the system will be at the ground state, i.e., 
without any phonons if there are no holes in the trap and 
in the polaronic state when there is one hole. Then tun- 
neling from the "L" lead will be tunneling to the state 
without any phonons. Therefore this process will probe 
only a single hole state, i.e. two step structure due to hole 
hopping. Tunneling to the "R" lead is from the ground 
state of the coupled hole-phonon system and so this tun- 
neling will produce the phonon steps in the TV depen- 
dence. Then if the gate voltage or the on-site energy 
is increased the tunneling to the "R" lead will become 
less sensitive to the internal structure of the hole-phonon 
state and hence the phonon steps would be suppressed. 
From Fig. 1 we conclude that for typical parameters of 
the DNA structure the hopping integral between the sites 
within the DNA traps and phonon frequency which de- 
termine the energetics of hole-phonon trap system, can 
be found from dependence of the tunneling current on 
Vsd- From a small gate voltage the phonon frequency 
can be found from the TV curve, while for a larger gate 
voltage the hopping integral can be obtained. 

The TV curve should show even richer structure for 
larger number of sites in the trap. In Fig. 2 the current 
as a function of bias voltage is shown for Nt = 3 sites. In 
this case tunneling is possible through the sites ia = I and 
«o = 2. For uncoupled hole-phonon system the TV curve 
shows three steps for iq — 1 (green hue in Fig. 2a), and 



two steps for ia — 2 (green line in Fig. 2b). This means 
that for iQ = 2 only two states have non-zero amplitude 
at i = 2 and contribute to the tunneling current. The 
finite hole-phonon coupling results in two effects: the 
phonon steps in the TV dependence similar to two-site 
trap (Fig. 1), and the polaronic effect which redistributes 
the hole density along the trap and increase or decrease 
the tunneling current. 

For small hole-phonon coupling (A — 0.5) the phonon 
steps are seen only at the first hopping step [Fig. 2a, b 
(black lines)]. The separation between the steps is the 
phonon frequency. Similar to Fig. 1, an increase of the 
gate voltage (on-site energy) suppresses the phonon steps 
and the TV curve becomes similar in structure to the 
uncoupled case [Fig. 2a, b (red Hues)]. For a larger hole- 
phonon coupling (A = 1) the steps due to hole hopping 
almost completely disappear for iq — I [Fig. 2c, black 
line], but some structure is still visible for io — 2 [Fig. 2d, 
black line] . For a larger gate voltage the steps due to the 
hopping overcome the phonon effects. 

The polaronic effects due to hole-phonon coupling can 
be clearly seen in Fig. 2. With incresing hole-phonon 
interaction the hole states become more localized at the 
center of the trap (see Fig. 3) which results in an increase 
of the current for tunneling through the central site of the 
trap io = 2 (see Fig. 2b, d), and decrease of the current 
through io = 1 (see Fig. 2a, c). In addition to changes to 
the tunneling current the polaronic effect also modifies 
the structure of the I-V curve. This can be seen in Fig. 2a, 
where with increasing hole-phonon interaction the third 
step due to the hole hopping disappears. 
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